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Methods based on the combination of the usual density functional theory (DFT) codes with 
the Hubbard models are widely used to investigate the properties of strongly correlated materi¬ 
als. Using first-principle calculations we study the electronic and magnetic properties of 20 half- 
metallic magnets performing self-consistent GGA+U calculations using both the atomic-limit (AL) 
and around-mean-field (AMF) functionals for the double counting term, used to subtract the corre¬ 
lation part from the DFT total energy, and compare these results to the usual generalized-gradient- 
approximation (GGA) calculations. Overall the use of AMF produces results similar to the GGA 
calculations. On the other hand the effect of AL is diversified depending on the studied material. In 
general the AL functional produces a stronger tendency towards magnetism leading in some cases 
to unphysical electronic and magnetic properties. Thus the choice of the adequate double-counting 
functional is crucial for the results obtained using the GGA+U method. 

PACS numbers: 75.50.Cc, 71.20.Lp, 71.15.Mb 


I. INTRODUCTION 

The rapid expansion of the field of spintronics and 
magneto-electronics brought magnetic materials at the 
nanoscale to the center of attention of modern electron¬ 
ics. The spin of the electron offers an additional degree of 
freedom in electronic devices with respect to conventional 
electronics based on semiconductors X The design of mag¬ 
netic nanomaterials with novel properties offers new func¬ 
tionalities to future devices, and to this respect ab-initio 
(also known as first-principles) studies of the electronic 
structure within density functional theory (DFT) play 
a crucial role allowing the modelling of the properties 
of several materials prior to their experimental growth. 
Among the most studied magnetic materials are the 
so-called half-metallic (HM) magnets*^ which present 
metallic behavior for the majority-spin electronic band 
structure and semiconducting for the minority-spin elec¬ 
tronic band structure. The ferromagnetic semi-Heusler 
compound NiMnSb was the first material for which the 
HM character was predicted and described*^ and since 
then several HM compounds have been discovered.— 
The implementation of half-metallic magnets in devices 
is an active field of research (see Ref. [8| for a review of 
the literature). 

DFT-based ab-initio electronic structure calcula¬ 
tions using either the local-spin-density approxima¬ 
tion (LSDA)^ or the generalized-gradient-approximation 
(GGA+i for the exchange-correlation functional are 
quite successful for magnetic materials from weak to in¬ 
termediate electronic correlations, but fail for systems 
with strong electronic correlations. There are two com¬ 
mon ways to include correlations in first-principles elec¬ 
tronic structure calculations. The first one is the so- 
called LDA+C7 scheme, in which the locaD(spin)-density 
approximation (L(S)DA) of DFT is augmented by an 
on-site Coulomb repulsion term and an exchange term 
with the Hubbard U and Hund exchange J parame¬ 


ters, respectively—^ Such a scheme has been applied 
for example to Co 2 FeSi, showing that correlations re¬ 
store the HM character of the compound,— and to 
NiMnSb— When the GGA functional is used instead 
of the L(S)DA the method is usually referred to as 
GGA+U scheme. A more elaborate modern computa¬ 
tional scheme, which combines many-body model Hamil¬ 
tonian methods with DFT, is the so-called LDA+DMFT 
method, where DMFT stands for Dynamical Mean-Field 
Theory j 15 : 16 LDA+DMFT has been applied to several 
HM magnetic systems like Co 2 MnSi*i£ NiMnSb+b— 
FeMnSb^i M^VAl^ VAs22 and CrAs : 24 : 25 

In the case of both LSDA+U (GGA+U) and 
LDA+DMFT schemes, the addition of the Hubbard U 
interaction introduces the need for a double-counting 
correction term in the energy functional to account for 
the fact that the Coulomb energy between the corre¬ 
lated states is already included in the LSDA (GGA) 
functional. Several double-counting schemes have been 
proposed in literature*^— and in all proposed schemes 
an averaged energy for the occupation of a selected ref¬ 
erence state is subtracted. Among the proposed func¬ 
tionals for the double-counting term, two are most com¬ 
monly used: the so-called around-mean-field (AMF) 
functional and the atomic-limit (AL) functional; the lat¬ 
ter is also referred to in literature as the fully local¬ 
ized limit (FLL) functional. The performance of these 
two functionals has attracted little attention in litera¬ 
ture. In 2009 Ylvisaker and collaborators presented an 
extensive study on the effect of the two functionals when 
performing self-consistent LSDA+U calculations for sev¬ 
eral magnetic materials— They have shown that the use 
of the LSDA+U interaction term usually enhances spin 
magnetic moments, but the AMF double-counting term 
gives magnetic states a significantly larger energy penalty 
than does the AL(FLL) functional and thus AL gives a 
stronger tendency to magnetism than AMF — 
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II. MOTIVATION AND COMPUTATIONAL 
METHOD 

As mentioned above, ab-initio electronic structure 
calculations based on the mixed LSDA+U/GGA+U 
schemes as well as LDA+DMFT are widely used to study 
the influence of electronic correlations on the electronic 
and magnetic properties of half-metallic magnets. Thus 
the study of the influence of the double-counting term on 
the calculated properties for these materials is extremely 
important with respect to their potential use in realistic 
devices. The aim of the present study is to explore the 
effect of both AL and AMF functionals when perform¬ 
ing GGA+U calculations with respect to usual electronic 
band structure calculations using the GGA functional 
for a wide range of half-metallic magnets (the reader is 
referred to Ref. Hg) for an extended discussion on the 
exact formulation of the two functionals). To achieve 
our goal we have employed the full-potential nonortlrog- 
onal local-orbital minimum-basis band structure scheme 
(FPLO) — For the GGA calculations we have used the 
Perdew-Burke-Ernzerhof parametrization— In the case 
of the GGA+U calculations the on-site Coulomb inter¬ 
actions for the correlated d or p orbitals are introduced 
via the Fo, F 2 , F 4 and Fq Slater parameters.— For all cal¬ 
culations a dense 20x20x20 grid in the reciprocal space 
has been used to carry out the integrals and both the 
charge density (up to 1CU 6 in arbitrary units) and the 
total energy (up to 10 -8 Hartree) have been converged 
in each case. 

In order to cover a wide range of half-metallic mag¬ 
nets in a coherent way, we have used in our calculations 
the ab-initio determined Coulomb effective interaction 
parameters (Hubbard U and Hund exchange J between 
localized d or p electrons) calculated in Ref. [H using the 
constrained Random Phase Approximation (cRPA) 33-36 
for 20 half-metallic magnets. We should note that (i) 
the determination of these parameters from experimen¬ 
tal data is a difficult task, and (ii) the constrained local- 
density approximation (cLDA), although is the most 
popular theoretical approach, 3 — — it is well known to give 
unreasonably large Hubbard U values for the late tran¬ 
sition metal atoms due to difficulties in compensating 
for the self-screening error of the localized electrons^ 
and thus cRPA which does not suffer from these dif¬ 
ficulties, although numerically much more demanding 
than cLDA, offers an efficient way to calculate the effec¬ 
tive Coulomb interaction parameters in solids. 33 ’- 38 We 
present results for all 20 half-metallic magnets studied 
in Ref. [32| which include representatives of the (i) semi- 
Heusler compounds like NiMnSb, (ii) ferrimagnetic full- 
Heusler compounds like Mn 2 VAl, (iii) inverse full-Heusler 
compounds like Cr 2 CoGa, (iv) usual L2i-type ferromag¬ 
netic full-Heusler compounds, (v) transition-metal pnic- 
tides like CrAs, and finally (vi) sp-electron (also called 
d°) ferromagnets like CaN. We have used the lattice pa¬ 
rameters presented in Table 1 of Ref. HD. The Slater 
parameters entering the FPLO method are connected to 


the Hubbard parameter Ulda+u and to the Hund ex¬ 
change J presented in Table II of Ref. [32| for the corre¬ 
lated p-states though the relations 

Fq = Ulda+u > F 2 = 5 x J, F4 = F 6 = 0, (1) 

and for the correlated d-states 

F 0 = U LD A + u , F *+ F * = ./,§■ = 0.625, F 6 = 0. (2) 

We should note here that Ulda+u is an effective param¬ 
eter depending both on the on-site intra-orbital Coulomb 
repulsion between electrons occupying the same orbital 
and on-site inter-orbital Coulomb repulsion between elec¬ 
trons occupying orbitals of the same £ character but dif¬ 
ferent mg value. Thus, our study covers a wide range 
of half-metallic magnets allowing for a deeper under¬ 
standing of the behavior of the AL and AMF double¬ 
counting functionals in the GGA+U calculated electronic 
and magnetic properties of different HM magnetic sys¬ 
tems. 


III. RESULTS AND DISCUSSION 
A. Binary Compounds 

We will start the presentation of our results from 
the binary compounds. There are two families of half- 
metallic binary compounds. The first includes the 
so-called sp-electron ferromagnets (also known as d°- 
ferromagnets) i 41 i 42 These compounds adopt the rocksalt 
cubic structure and have no transition-metal atoms in 
their chemical formula. We consider the nitrides and the 
carbides (CaN, SrN, SrC, and BaC) since they have the 
largest calculated Curie temperatures among the stud¬ 
ied sp-electron ferromagnets.— Their total spin mag¬ 
netic moment in units of ps equals 8 — Z t , where Z t is 
the total number of valence electrons in the unit cell; 
this behavior is known as Slater-Pauling (SP) — A de¬ 
tailed discussion on the origin of this rule and its con¬ 
nection to the half-metallicity can be found in Ref. [42]. 
The usual GGA calculations produced for all four stud¬ 
ied compounds a half-metallic state with total spin mag¬ 
netic moments of 1 /is for the nitrides and 2 ps for the 
carbides. The results are gathered in Table Q] The spin 
moment is carried mainly by the N and C atoms. Our cal¬ 
culated GGA results are similar to the GGA ones derived 
with the full-potential linearized augmented plane-wave 
(FLAPW) method as implemented in the FLEUR 3 - code 
in Ref. [32j. The use of the AMF within the GGA+U 
scheme leaves intact both the calculated spin magnetic 
moments and density of states (DOS) with respect to 
GGA calculations (we do not present the DOS since they 
are similar to the ones presented in literature). On the 
contrary the use of the AL functional has a tremendous 
effect on the calculated results. It produces an unreason¬ 
able and unphysical charge transfer from the Ca(Sr,Ba) 
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TABLE I: Atom-resolved and total spin magnetic moment 
per formula unit for the XY binary compounds. Results have 
been obtained within the FPLO method^ using the GGA 
functional for the exchange interaction potential and the 
GGA+U scheme employing both the atomic-limit (AL - also 
known as fully-localized-limit FLL) and the around-mean- 
field (AMF) functionals for the double counting term. Val¬ 
ues for the on-site Coulomb and exchange parameters are the 
ab-initio determined ones within the constrained Random- 
Phase-Approximation (cRPA) in Ref. [32|. Lattice constants 
are the ones presented in Table 1 in the later reference. Note 
that for the compounds which do not contain transition metal 
atoms (known as d°-ferromagnets) GGA+U within the AL 
functional gives unrealistic results. 


Comp. 

Functional 

m x 

rriY 

fT^total 

CaN 

GGA 

-0.065 

1.065 

1.000 


GGA+U (AL) 

11.836 

-4.836 

7.000 


GGA+U (AMF) 

-0.065 

1.065 

1.000 

SrN 

GGA 

-0.072 

1.072 

0.999 


GGA+U (AL) 

16.363 

-9.362 

7.000 


GGA+U (AMF) 

-0.072 

1.072 

0.999 

SrC 

GGA 

-0.004 

2.004 

1.999 


GGA+U (AL) 

15.578 

-9.578 

6.001 


GGA+U (AMF) 

-0.004 

2.004 

1.999 

BaC 

GGA 

0.057 

1.943 

2.000 


GGA+U (AL) 

3.261 

0.378 

3.999 


GGA+U (AMF) 

0.057 

1.943 

2.000 

VAs 

GGA 

2.427 

-0.427 

2.000 


GGA+U (AL) 

2.415 

-0.415 

2.000 


GGA+U (AMF) 

2.151 

-0.151 

2.000 

CrAs 

GGA 

3.614 

-0.614 

3.000 


GGA+U (AL) 

3.880 

-0.880 

3.000 


GGA+U (AMF) 

3.541 

-0.541 

3.000 

MnAs 

GGA 

4.173 

-0.311 

3.862 


GGA+U (AL) 

4.476 

-0.476 

3.999 


GGA+U (AMF) 

3.979 

-0.326 

3.652 


atoms to the N(C) atoms resulting to huge values of the 
atom-resolved spin moments. This state is obviously an 
artifact of the method. We cannot explain the origin of 
this behavior but starting form various configurations all 
calculations involving the AL functional converged to the 
same results and thus the breakdown of the AL should 
be attributed to its characteristics. 

The second family of binary compounds under study 
are the binary VAs, CrAs, and MnAs transition metal 
pnictides. The first observation of such a compounds 
being half-metal was made in 2000 when Akinaga 
and his collaborators managed to grow multilayers of 
CrAs/GaAs^ CrAs was found to adopt the zincblende 
structure of GaAs and was predicted to be a half-metal 
with a total spin magnetic moment of 3 hb in agreement 
with experiments^ Several studies followed this initial 
discovery, and electronic structure calculations have con¬ 
firmed that also similar binary XY compounds, where X 
is an early transition-metal atom and Y an sp element, 
should be half-metals and the total spin magnetic mo¬ 
ment follows a SP rule similar to d-ferromagnets being 
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FIG. 1: (color online) Total density of states (DOS) as a 
function of the energy for the CrAs compound within the 
GGA+U method using both the atomic-limit (AL) and the 
around-mean-field (AMF) functionals for the double-counting 
term. GGA+U results are compared to the GGA calculated 
DOS. The zero in the energy axis has been set to the Fermi 
level. Positive(negative) values of the DOS correspond to the 
majority(minority)-spin electrons. 


now equal to Z t - 8 M& 

In Table [T| we gathered all the calculated spin mag¬ 
netic moments. GGA gives a half-metallic state for VAs 
and CrAs, while for MnAs the Fermi level is slightly 
above the minority-spin energy gap and the total spin 
magnetic moment slightly smaller than the ideal value 
of 4 pb for half-metallicity to occur. These results have 
been largely discussed in literature . 54 ! 55 For both VAs 
and CrAs, GGA+U self-consistent calculations yield a 
half-metallic state within both AL and AMF function¬ 
als with the same total spin magnetic moment but with 
substantial variations of the atom-resolved spin magnetic 
moments. For MnAs the use of AL functional leads to a 
half-metallic state contrary to AMF for which the Fermi 
level is above the minority-spin gap. Overall AL leads to 
larger absolute values of the atomic spin moments with 
respect to GGA while AMF leads to smaller values. This 
behavior of the atomic spin magnetic moments confirms 
the conclusion in Ref. |29| that AMF gives the magnetic 
state a large energy penalty with respect to AL. 

Since DOS present similar trends between the three 
transition metal binary compounds, we present in Fig. Q] 
the calculated DOS per formula unit for CrAs. In the 
upper panel we compare the GGA+U calculated DOS 
within the AL functional to the usual GGA calculated 
DOS, and in the lower panel we present a similar graph 
for the AMF case. In the presented energy Cr DOS dom- 





































inates. GGA produces a large minority-spin gap with a 
large exchange splitting between the occupied majority- 
spin bands and the unoccupied minority-spin bands and 
thus strong tendency to magnetism manifested also by 
the large (~3.6 hb) Cr spin moment. The use of the AL 
double-counting functional in the GGA+U calculations 
lead to an almost rigid shift of the minority spin DOS 
towards higher energies, while in the majority-spin DOS 
only the double-degenerate e g states at about -1.5 eV 
move lower in energy (see Ref. [54] for a discussion of the 
character of the bands). In the case of AMF the majority- 
spin band structure shows a similar behavior with respect 
to the GGA results as the AL case. But in the minority- 
spin band structure the tendency is the opposite now. 
Since AMF does not favor magnetism as strongly as AL, 
the minority-spin band structure now presents an almost 
rigid shift towards lower energy values. These finding 
also explain the behavior of the MnAs compound. In 
the case of AL the minority-spin band structure moves 
towards higher energy values and the Fermi level now 
moves within the gap and half-metallicity appears. 


B. Heusler compounds 

Heusler compounds are a huge family of intermetallic 
compounds presenting various types of electronic and 
magnetic behaviors ^-* — Several among them are half- 
metallic ferromagnets / ferrimagnets / antiferromagnets 
and are of particular interest due to their very high 
Curie temperatures, which usually exceed 1000 K, mak¬ 
ing them ideal for applications^ There are four main 
families of Heusler compounds: (i) the semi-Heuslers 
also known as half-Heuslers like NiMnSb which have the 
chemical type XYZ with X and Y being transition metal 
atoms, (ii) the usual full Heuslers like Co 2 MnSi with the 
chemical type X 2 YZ where the valence of X is larger than 
the valence of Y and the two X atoms are equivalent, 
(iii) the quaternary Heuslers like (CoFe)MnSi which 
present similar properties with the full-Heuslers, and 
finally (iv) the so-called inverse-Heuslers, like C^CoGa 
which have also the chemical type X 2 YZ but now the 
valence of X is smaller than the valence of Y and due to 
the change of the sequence of atoms in the unit cell the 
two X atoms are no more equivalent^ i§Z We present 
results for all families of compounds with the exception 
of quaternary-Heuslers which present similar behavior to 
the full-Heuslers and for which no Hubbard parameters 
have been derived in Ref. [32). 


1. Semi-Heuslers 

The first family of Heusler compounds for which we 
will present results are the semi-Heuslers. The first com¬ 
pound that was predicted to be a half-metal was actu¬ 
ally a semi-Heusler, NiMnSb^ Their total spin magnetic 
moment follows also a SP rule being Z t -18 (for an ex¬ 



FIG. 2: (color online) Ni and Mn atom-resolved DOS in 
NiMnSb. Details as in Fig. |Tj 


tended discussion see Ref. [5§). In Table HD we have gath¬ 
ered our calculated spin magnetic for all studied cases 
and for three compounds FeMnSb, CoMnSb and NiMnSb 
(note that for FeMnSb we were not able to converge the 
GGA+U calculations using the AMF functional). As we 
move from one compound to the other, the total num¬ 
ber of valence electrons increases by one and so does the 
GGA calculated total spin magnetic moment. Mn atoms 
in all case posses a large value of spin magnetic moment 
which starts from ~3.4 /is in FeMnSb and exceeds 4 /is 
in NiMnSb. As we increase the total number of valence 
electrons the spin magnetic of the X atoms also increases 
being ~-1.3 /zb for Fe, -0.34 /zb for Co and 0.14 /zb for 
Ni in the corresponding compounds. The GGA calcu¬ 
lated DOS, presented in Fig. [2] for NiMnSb, has been 
studied in detail in literature and it is mainly character¬ 
ized by the large exchange splitting between the occupied 
majority-spin and the unoccupied minority-spin d-states 
at the Mn site which together with the very small weight 
of the occupied minority-spin states are responsible for 
the large Mn spin magnetic moments. This feature is 
common for all three studied compounds and has been 
already observed in literature.— 

The self-consistent GGA+U calculations using the 
AMF functional for the CoMnSb and NiMnSb com¬ 
pounds produced a similar picture to the GGA calcu¬ 
lations. The total spin magnetic moment, as shown in 
Table ITT1 remains identical to the GGA case and the atom- 
resolved spin magnetic moments only scarcely changed. 
This is also reflected on the Ni and Mn resolved DOS for 
NiMnSb in Fig. 0 where the GGA+U within AMF cal¬ 
culated DOS is almost identical to the GGA calculated 
DOS. The effect of the use of the AL functional is more 
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TABLE II: Similar to Table[I]for the semi-Heusler compounds 
crystallizing in the Clt lattice structure having the XYZ 
chemical formula. 


Comp. 

Functional 

m x 

rriY 

mz 

^total 

FeMnSb 

GGA 

-1.279 

3.374 

-0.095 

2.000 


GGA+U(AL) 

-2.274 

4.559 

0.042 

2.327 

CoMnSb 

GGA 

-0.340 

3.568 

-0.227 

3.000 


GGA+U (AL) 

-1.264 

4.520 

-0.186 

3.068 


GGA+U (AMF) 

-0.358 

3.535 

-0.177 

3.000 

NiMnSb 

GGA 

0.143 

4.031 

-0.174 

3.999 


GGA+U (AL) 

-0.087 

4.553 

-0.301 

4.164 


GGA+U (AMF) 

0.144 

4.026 

-0.171 

3.999 


drastic. As shown in Fig. [2] GGA+U within the AL 
functional compared to usual GGA leads to large mod¬ 
ifications of the DOS of the transition metal atoms. In 
the case of Mn, the exchange splitting between the oc¬ 
cupied majority-spin and the unoccupied minority-spin 
states increase considerably. In the majority-spin band 
structure of Mn, the occupied states move lower in en¬ 
ergy and as a result they are no more in the same en¬ 
ergy with the Ni majority-spin states. This leads to a 
weaker hybridization between the d states of Ni and Mn 
atoms and in the Ni DOS the width of the majority bands 
becomes smaller as a result of the weaker hybridization 
effects. Almost all the weight of the occupied minority- 
spin band structure is located at the Ni atoms while al¬ 
most all unoccupied minority-spin states are located at 
the Mn atom. Thus there is almost no hybridization be¬ 
tween the minority-spin d-states of Ni and Mn and the 
former are not affected by the shift of the later to larger 
energy values being identical to the GGA case. These 
changes in DOS are also reflected on the spin magnetic 
moments in Table EH The larger tendency to magnetism 
within AL compared to AMF leads to slightly larger to¬ 
tal spin magnetic moments which deviate from the ideal 
integer values of the SP rule and the Fermi level is lo¬ 
cated slightly below the minority-spin energy gap. In the 
case of FeMnSb and CoMnSb, both the absolute values 
of the Fe(Co) and Mn spin magnetic moments increase 
by about 1 hb almost cancelling each other. In the case 
of NiMnSb the variations in the atomic spin magnetic 
moments are considerably smaller since almost all Ni d- 
states are occupied in all studied cases. But even for 
NiMnSb the Mn spin moment increases by ~0.5 hb and 
the Ni spin moment decrease by about 0.3 fi b■ The Sb 
atoms in all cases also present changes in their atomic 
spin magnetic moments between the various calculations 
although these variations are considerably smaller than 
for the transition metal atoms. 


2. Full-Heuslers 

The second family of Heusler compounds which may 
present half-metallicity are the so-called usual full- 
Heuslers crystallizing in the cubic L2i structures. Half- 


TABLE III: Similar to Table|]]for the full-Heusler compounds 
crystallizing in the L2i lattice structure having the X 2 YZ 
chemical formula. 


Comp. 

Functional 

m x 

rriY 

mz 

^total 

Mn 2 VAl 

GGA 

-1.670 

1.227 

0.113 

-1.999 


GGA+U (AL) 

-4.102 

3.169 

0.483 

-4.551 


GGA+U (AMF) 

-1.798 

1.527 

0.079 

-1.990 

Mn 2 VSi 

GGA 

-0.801 

0.557 

0.060 

-0.985 


GGA+U (AL) 

-2.248 

2.759 

0.376 

-1.362 

Co 2 CrAl 

GGA 

0.737 

1.684 

-0.160 

2.999 


GGA+U (AMF) 

0.965 

1.222 

-0.153 

3.000 

Co 2 CrSi 

GGA 

0.934 

2.242 

-0.111 

4.000 


GGA+U (AL) 

0.871 

2.525 

-0.267 

4.000 


GGA+U (AMF) 

0.890 

2.187 

0.031 

4.000 

Co 2 MnAl 

GGA 

0.673 

2.910 

-0.231 

4.025 


GGA+U (AL) 

1.381 

4.176 

-0.501 

6.438 


GGA+U (AMF) 

1.048 

1.998 

-0.096 

3.999 

Co 2 MnSi 

GGA 

0.972 

3.195 

-0.140 

4.999 


GGA+U (AL) 

0.732 

3.954 

-0.338 

5.080 


GGA+U (AMF) 

0.987 

3.020 

-0.004 

5.000 

Co 2 FeAl 

GGA 

1.163 

2.870 

-0.203 

4.999 


GGA+U (AL) 

1.188 

3.326 

-0.520 

5.177 


GGA+U (AMF) 

1.216 

2.673 

-0.105 

4.999 

Co 2 FeSi 

GGA 

1.327 

2.926 

-0.042 

5.539 


GGA+U (AL) 

1.375 

3.450 

-0.201 

5.999 


metallicity can be combined either with the appearance 
of ferrimagnetism, when the X atoms in the X 2 YZ is the 
Mn one, or with ferromagnetism when X is Co. In all 
cases the total spin magnetic moment in /is follows a 
SP rule being Z t - 24i£i In Table HTT1 we have gathered the 
calculated spin magnetic moments for all studied com¬ 
pounds with both GGA and GGA+U methods using 
both the AL and AMF double-counting functional in the 
later case. When one case is missing in the table, this is 
due to the fact that we were not able to get convergence 
irrespectively of the starting input which we have used. 

First, we will discuss our results on the half-metallic 
ferrimagnetic Mn 2 VAl and Mn 2 VSI compounds where 
the total spin magnetic moments is negative since the 
total number of valence electrons is less than 24. More¬ 
over the Mn spin magnetic moments are antiferromag- 
netically coupled to the V spin moments due to their 
small distanced In the case of Mn 2 VAl, GGA+U cal¬ 
culations within AMF produced similar spin moments 
and DOS to the GGA case; we were not able to con¬ 
verge GGA+U within AMF for the Mn 2 VSi compound. 
For both compounds the use of AL double-counting func¬ 
tional produced unpliysical results similar to the case of 
d°-ferromagnets in Sec. Emu The use of AL tripled, 
with respect to the GGA case, the absolute values of the 
spin magnetic moments of the transition metal atoms in 
Mn 2 VA; in the case of V in Mn 2 VSi the increase is almost 
600% . Thus the use of AL for the half-metallic ferrimag¬ 
netic Heusler compounds obviously is inadequate. 

In the case of ferromagnetic full-Heuslers containing 
Co the effect of using both AMF and AL on the calcu¬ 
lated electronic and magnetic properties is more complex 























FIG. 3 : (color online) Co and Mn atom-resolved DOS in 
C02M11SL Details as in Fig. [l] 


than in all the previously studied cases. When Y is Cr 
(Co 2 CrAl and Co 2 CrSi) both AL and AMF yielded a per¬ 
fect half-metallic state with the total spin magnetic mo¬ 
ment being equal to the ideal values predicted by the SP 
rule as shown in Table Hill When Y is Mn (Co 2 MnAl and 
Co 2 MnSi) AMF produced a half-metallic states and both 
atom-resolved and total spin magnetic moments where 
close to the GGA case, but AL led to a considerable 
increase of the Mn spin moment similarly to the semi- 
Heuslers. The increase of the Mn spin moment within 
AL led to an increase also of the total spin magnetic mo¬ 
ment which is only 0.80 ^lb for Co 2 MnSi but reaches 
the ~1.4 hb for Co 2 MnAl. When Y is Fe (Co 2 FeAl 
and Co 2 FeSi) the behavior of the spin moments with re¬ 
spect to the GGA results is similar within both AL and 
AMF to the case where Y is Mn. Moreover in the case 
of Co 2 FeSi which is not half-metallic within GGA, the 
use of GGA+U combined with AL leads to a total spin 
magnetic moment of 6 ur and to a half-metallic state 
as shown also in Ref. [l3|. Although the GW scheme^ 
produced similar results to the GGA+ calculations, cor¬ 
relations in this materials are still an open issue since 
recent results by Meinert and collaborators show that a 
self-consistent calculation fixing the total spin magnetic 
moment to 6 reproduces more accurately the posi¬ 
tion of the band with respect to available experimental 
dataM 

To understand the behavior of the spin moments we 
have to examine in detail the behavior of the DOS. Since 
the trends when either AMF or AL is employed are simi¬ 
lar for all six ferromagnetic Co-based full-Heuslers under 
study, we will use Co 2 MnSi as an example and in Fig. 
[3] we present the Co and Mn resolved DOS. When the 



FIG. 4 : (color online) Total DOS per formula unit for the 
Co2MnSi and Co2MnAl compounds. Details as in Fig. [T] 


GGA+U combined with AMF is used (left panel) there is 
a significant change in the DOS unlikely all other families 
of half-metallic compounds discussed above. AMF en¬ 
hances the tendency to magnetism with respect to GGA. 
For the Mn atom the occupied majority spin states shift 
lower in energy but the minority-spin energy gap in the 
Mn DOS remains unchanged. Through hybridization 
also the Co majority-spin DOS shifts lower in energy 
and so do also the occupied Co minority-spin states. This 
leads to an increase of the energy gap in the Co minority- 
spin band structure. Since as explained in Ref. [6l| Co 
atoms present in usual GGA calculations a much smaller 
gap than the Mn atoms and thus determine the energy 
gap in the total DOS, the opening of the former further 
stabilizes the half-metallic state. 

The GGA+U method combined with AL even fur¬ 
ther enhances the tendency to magnetism with respect to 
AMF as concluded in Ref. [29l In the case of Mn atoms 
the exchange splitting between the occupied majority- 
spin and unoccupied minority-spin states is greatly en¬ 
hanced as for Mn in NiMnSb and thus the energy gap 
becomes much larger. As a side-effect some weight in 
the minority-spin band structure appears just below the 
Fermi level. Thus although with respect to the GGA 
case, AL opens the gap the Fermi level is located close to 
the left edge of the gap instead of the middle in the GGA 
case. Co DOS follows through hybridization the behavior 
of the Mn d-states and the gap is now also much larger 
but the occupied minority-spin states move closer to the 
Fermi level which now just crosses the states just below 
the low-energy edge of the gap and the total spin mag¬ 
netic moment within AL is slightly larger than the ideal 
value of 5 n b ■ 
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TABLE IV: Similar to Table |I] for the full-Heusler compounds 
crystallizing in the inverse XA lattice structure having the 
Cr 2 YZ chemical formula, where the two Cr atoms occupy sites 
of different symmetry (see text). 


Comp. 

Functional m Cr A 

m Cr B 

my 

mz 

^total 

Cr 2 FeGe 

GGA -1.454 

1.753 

-0.313 

0.027 

0.012 


GGA+U(AL) -4.162 

4.672 

-1.054 

0.551 

0.006 

Cr 2 CoGa GGA -2.680 

1.973 

0.379 

-0.014 

0.069 


GGA+U (AL) -4.860 

4.137 

1.160 

-0.082 

0.520 


In the case of Co 2 MnAl the change in the spin mag¬ 
netic moments is larger within both AL and AMF func¬ 
tionals. As shown in Fig. [2 although we just change 
A1 for Si in Co 2 MnSi, the AMF DOS shows a different 
tendency with respect to the energy gap. The exchange 
splitting between occupied majority and unoccupied mi¬ 
nority spin states is smaller, and within AMF the gap 
is smaller than within GGA showing the contrary ten¬ 
dency to Co 2 MnSi where AMF produced a larger gap 
with respect to GGA. For Co 2 MnAl within usual GGA 
the Fermi level is close to the left edge of the gap while 
for Co 2 MnSi it is located at the middle of the gap. Thus 
in the case of AL based calculations the shift of the Co 
occupied minority spin states towards higher energies for 
Co 2 MnAl, discussed just above also for Co 2 MnSi, leads 
to the loss of the half-metallicity since now the Fermi level 
crosses the occupied minority-spin states. The other Al- 
based Heuslers (Co 2 CrAl and Co 2 FeAl) exhibit within 
GGA a DOS around the minority-spin energy gap simi¬ 
lar to Co 2 MnSi and not Co 2 MnAl and thus the increase 
in their total spin magnetic moment within AL is much 
smaller than for Co 2 MnAl. 

3. Inverse-Heuslers 

The last family of potential half-metallic Heusler com¬ 
pounds ar the so-called inverse Heusler compounds^ 
Among these half-metals the most interesting are the 
so-called fully-compensated ferrimagnets (also known 
as half-metallic antiferromagnets) like C^FeGe and 
Cr 2 CoGa. These materials are of special inter¬ 
est since they combine half-metallicity to a zero to¬ 
tal net magnetization and thus are ideal for spin- 
tronic/magnetoelectronic devices due to the vanishing 
external stray fields created by them^ We should note 
that films of C^CoGa have been grown experimentally^ 
and this compounds has been predicted to exhibit ex¬ 
tremely large Curie temperature^ 7 As shown in Table [IV] 
GGA yields for both Cr 2 FeGe and Cr 2 CoGa compounds 
a total spin magnetic moment close to zero (for an ex¬ 
tended discussion on the half-metallic inverse Heuslers 


see Ref. [(34). Note that we have two inequivalent Cr 
atoms in these compounds denoted by the superscripts 
A and B in Table [TV] We were not able to converge 
the GGA+U self-consistent calculations using the AMF 
double-counting functional. For the AL functional al¬ 
though the total spin magnetic moment stays close to 
zero, the absolute values of the Cr spin magnetic mo¬ 
ments are about doubled leading to an unphysical sit¬ 
uations. Thus for these materials the use of GGA+U 
combined with AL is not able to produce a reasonable 
description of the electronic structure as was also the case 
for the semi-Heuslers and the ferrimagnetic full-Heuslers. 


IV. CONCLUSIONS 

We have studied the electronic and magnetic proper¬ 
ties of 20 half-metallic magnets performing self-consistent 
GGA+U calculations using both the atomic-limit (AL) 
and around-mean-field (AMF) functionals for the dou¬ 
ble counting term and compared them to the usual GGA 
calculations. Overall the use of AMF produced results 
similar to the usual GGA calculations. The effect of AL 
was diversified depending on the studied material. In 
the case of d°-ferromagnets, semi-Heuslers, ferrimagnetic 
full-Heuslers and inverse Heuslers the use of AL leads to 
unrealistic electronic and magnetic properties of the stud¬ 
ied compounds and thus its use is not justified. On the 
other hand in the case of transition-metal binary com¬ 
pounds and usual ferromagnetic full-Heusler compounds 
the use of AL enhanced the tendency towards magnetism 
with respect to both GGA and GGA+U combined with 
AMF. Depending on the position of the Fermi level, there 
were cases like MnAs and Co 2 FeSi for which AL produced 
a half-metallic state contrary to GGA and GGA+U com¬ 
bined with AMF, cases like VAs, CrAs and Co 2 CrSi 
where all three methods produced a half-metallic state, 
and cases like Co 2 MnAl, Co 2 MnSi and Co 2 FeAl where 
the use of AL led to the loss of half-metallicity. 

Methods based on the combination of the usual density 
functional theory (DFT)-based codes and of the Hubbard 
U- Hund’s exchange J are widely used to investigate the 
properties of strongly correlated materials. Our results 
suggest that especially in the case of half-metallic mag¬ 
nets the choice for the double counting functional used 
to subtract the part from the DFT total energy, which 
is associated to the Coulomb repulsion between the cor¬ 
related orbitals, plays a decisive role on the obtained re¬ 
sults. Thus the study of correlations in half-metallic mag¬ 
nets is still an open issue and further studies are needed 
to establish the predictive power of methods based on 
the U and J parameters like GGA+U or LDA+DMFT 
methods. 
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